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We track the secondary bifurcations of coherent states in plane Couette flow and show that they 
undergo an incomplete periodic doubling cascade that ends with a crisis bifurcation. We introduce 
a symbolic dynamics for the orbits and show that the ones that exist fall into the universal sequence 
described by Metropolis, Stein and Stein for unimodal maps. The periodic orbits cover much of the 
turbulent dynamics in that their temporal evolution overlaps with turbulent motions when projected 
onto a plane spanned by energy production and dissipation. 



The identification of the different 'routes to 
chaos' [H [2] has suggested ways in which they 
can also be 'routes to turbulence' that lead from 
a laminar flow state through various bifurcations 
to a turbulent attractor. Studies of the turbu- 
lence transition in plane Couette flow and pipe 
flow, which do not have a linear instability of 
the laminar profile, have revealed a route to tur- 
bulence that passes through the formation of a 
chaotic saddle [SlS]. The saddle develops around 
a scaffold of coherent structures that appear in 
saddle-node type bifurcations, with the notable 
difference to the standard scenario that also the 
'node'-states are typically unstable. In this tran- 
sitional region the turbulence is not persistent 
but transient, and can hence be associated with a 
transient chaotic saddle rather than an attractor. 
In the spirit of a dynamical systems approach to 
characterizing these features [5j, we here present 
an analysis of the phase space structures near 
transition in plane Couette flow and study the 
various bifurcations and periodic orbits that ap- 
pear. 

I. INTRODUCTION 

The application of dynamical systems ideas to fluid 
mechanical problems has provided much insight into the 
ways fluids transit from a laminar state to a turbulent 
one. The various bifurcations that appear in the case 
of fluids heated from below (Rayleigh-Benard) or fluids 
driven by centrifugal instabilities (Taylor- Couette) have 
been studied experimentally and theoretically in consid- 
erable detail [6j [7] . Studies of turbulence in liquid he- 
lium have shown the period doubling cascade [8] and the 
quasi-periodic routes to chaos [9 , and the rich structures 
of secondary bifurcations and flow states in RB and TC 
have been explored [6] [7J [10] . The approach has been so 
successful that flows without linear stability of the linear 
profile, like plane Couette flow and pipe flow, have been 



suspected to belong to the class of subcritical bifurca- 
tions [11 , perhaps with the bifurcation point at infinity 
(as reflected, e.g., in the title of Nagatas first papers on 
coherent structures in pipe flow [12]). However, the first 
structures that are observed near the transition point are 
not steady or simple flows, but they are fluctuating spa- 
tially and temporally, indicating that these states un- 
dergo many further bifurcations as the Reynolds number 
is lowered. 

A more promising approach starts with the turning 
points of the subcritical bifurcations, the saddle-node bi- 
furcations in which these structures appear. The possi- 
bility to approach flows from this point of view is sug- 
gested by the remarkable persistence of key features of 
these structures [3 : they are all dominated by down- 
stream vortices which modify the base profile to gener- 
ate streaks. This by itself is not enough, as flows that 
are translationally invariant in the downstream direction 
cannot be sustained [13]. So three-dimensionality is an- 
other key element, and also the reason for the absence of 
analytical studies of this transition scenario. 

The importance of vortices and streaks has been de- 
duced from numerical studies of plane Couette flow in 
small domains, where a recurrent pattern of vortices and 
streaks could be identified [14] [15]. The interactions 
between vortices and streaks and the instabilities that 
streaks can undergo has been studied by Waleffe [16], 
who proposed a self-sustaining cycle (analyzed also in 
[T7H2Q] ). The cycle is usually explained by starting with 
the vortices, which then drive streaks, and the streaks 
then undergo an instability in which normal vorticity is 
created. To close the cycle, this normal vorticity is then 
argued to be turned into downstream vorticity by the 
mean profile. 

An example where the transitions between the differ- 
ent states and their symmetry can he analyzed directly, 
and in a physically meaningful setting, is Taylor- Couette 
flow in the limit in which it approaches plane Couette 
flow [21 . In the curved case, say with the inner cylinder 
rotating and the outer one at rest, the flow undergoes a 
first bifurcation to Taylor-vortex flow, in which vortices 
with an azimuthal symmetry appear. As the rotation 
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rate is increased, the vortices undergo a secondary bi- 
furcation, in which an azimuthal modulation is added 
and the rotational symmetry is broken. Increasing the 
radii while keeping the distance between the cylinders 
fixed, the influence of the curvature becomes smaller and 
the flow approaches plane Couette flow. Along the way, 
the bifurcation structure changes: the transition point to 
the formation of azimuthally symmetric vortices moves 
out infinity, so that the linear stability of plane Couette 
flow is recovered. The three-dimensional wave structures, 
however, can be traced over to plane Couette flow, where 
they now appear in a saddle node bifurcation [21] . In go- 
ing from the Taylor- Couette to plane-Couette flow the 
transition can be studied within realistic models. In the 
case of pipe flow no physical family of flows could be iden- 
tified, so a body force was introduced to embed the sys- 
tem into one where the bifurcations can be realized. De- 
spite this artificial setting, the coherent structures thus 
found are similar to the ones in plane Couette flow and 
share many properties [22] . 

The coherent structures just described are usually un- 
stable, and can appear only transiently in the flow [23]- 
[25] . Moreover, it was found that they are not sufficiently 
entangled to close off the phase space around them to 
form an attract or: thus it is not only the coherent struc- 
tures that are visited transiently, but the entire turbulent 
dynamics becomes a transient phenomenon, character- 
ized by an exponential distribution of lifetimes, as can 
be expected for a chaotic saddle [26, 27 J. While this sets 
the stage from the dynamical systems point of view, a 
couple of questions remain to be addressed [4]. We will 
address two of them, related to the secondary bifurca- 
tions of coherent structures, in the following. 

It was noticed long ago in [28] , that the coherent struc- 
tures that appear in plane Couette flow at Re= 125 are 
stable very close to the bifurcation point and become un- 
stable via Hopf bifurcations. This would suggest that the 
transition creates an attractor rather than a transient 
saddle. However, at higher Re the flowed evidently is 
transient: so how does the phase space near these coher- 
ent states change as one increases the Reynolds number 
and how does the flow change from chaotic but persistent 
to transient? The analysis of these properties in pipe flow 
shows that the sequence of bifurcations that lead from a 
(small) chaotic attractor to transient dynamics can be 
rather intricate [29] [30] . 

A second item concerns the possibility to describe the 
flows using periodic orbits, as is possible in other dy- 
namical systems [5] [31] [32] . Identifying periodic orbits in 
high-dimensional systems is a major technical challenge, 
which has been overcome in a few cases only. Most no- 
tably, for the cases studied here, a few isolated periodic 
orbits have been found in [33 . In the present case their 
identification is considerably simplified by the presence 
of a period doubling cascade. 

The outline of the paper is as follows: in section [TT] 
we discuss the organization of the phase space near the 
lowest saddle-node bifurcation in plane Couette flow. In 



section Hill we discuss the characterization of the internal 
dynamics and bifurcations in terms of a discrete map. We 



follow this in section IV with a discussion of the periodic 
orbits and their properties. 



II. BIFURCATIONS 

We study an incompressible Newtonian fluid between 
two infinitely extended plates in the x-z plane located at 
y = ±h that move with speed ±Uo along the x-axis. A 
length-, velocity and time-scale are given by h, Uq and 
t — h/Uo. The control parameter is the Reynolds number 
Re = Uod/v, where v is the viscosity of the fluid. We use 
periodic boundary conditions in the spatially extended 
directions and no-slip boundary conditions at the walls, 
with a box size of 2tt x 2 x tt in the downstream, wall- 
normal and spanwise directions, respectively. We use 
channelflow [34] to integrate this system with a Fourier- 
Chebyshev-Fourier scheme and a numerical resolution of 
32 x 33 x 32 modes. The main quantity we use in this 
study is the cross-flow energy, i.e. the amount of energy 
in the flow components transverse to the laminar base 
flow, 
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or the square root of E c f so as to reduce the range of it. 

The box size in this study is just slightly longer and 
narrower than the cell studied in Ref. [15 . In this and 
other smaller boxes the typical flow consists of a pair of 
counter-rotating streamwise vortices and a pair of alter- 
nating streaks. This limitation in the number of possible 
flow structures greatly simplifies the sequence of bifurca- 
tions. To simplify the analysis even further we enforce a 
shift- and-reflect symmetry 

a[u,v,w](x,y,z) = [u, v, -w](x + L x /2, y, -z), (2) 

which, up to a discrete displacement by Lz/2, fixes the 
spanwise location of the flow structures. 

The different box size compared to the calculations by 
Nagata [12] and Clever and Busse [28] has the effect of 
moving the first bifurcation to a higher Reynolds number, 
namely Re = 163.8, where a pair of fixed points appears 
in a saddle-node bifurcation. The lower branch equilib- 
rium (LB) has one unstable direction, the upper branch 
equilibrium (UB) is linearly stable. Since the laminar 
profile is linearly stable for all Reynolds numbers [35] , 
there now coexist two attractors. The separatrix of their 
basins of attraction is formed by the stable manifold of 
LB. 

The situation is illustrated in figure [TJ a). The initial 
conditions that correspond to the points in the plane of 
the figures are spanned by the laminar profile, and the 
upper and lower branch states. The x-axis is a linear in- 
terpolation between LB and UB and the y-axis gives the 



3 



a) Re = 164 




3 

a 

O 

CD 

PU 



b) Re = 188 



0.42 I 

0.40 

0.38 

0.36 

0.34 

0.32 

0.30 

0.28 

0.26 



-0.4 



o 



-0.2 0.0 



A 



o 



0.2 



0.4 



0.6 



c) Re = 195 



0.42 









a 



1000 

900 

800 

700 

600 

500 

400 

300 

200 

100 



1000 

900 

800 

700 

600 

500 

400 

300 

200 

100 



1000 

900 

800 

700 

600 

500 

400 

300 

200 

100 



perturbation amplitude with which the state is scaled. 
Formally, the initial conditions are given by 



FIG. 1. The basins of attraction of the attractors in plane 
Couette flow for three Reynolds numbers near the bifurcation 
point of the Nagata-Clever-Busse-state. Basins are depicted 
by the lifetime of initial conditions. States on the x-axis are 
obtained by interpolating between the lower- and the upper- 
branch solution, rescaled by a factor that ensures that the 
distance between the two solutions in the plot is equal to the 
L 2 -norm of their difference, see equation Q. These states are 
then scaled in amplitude along the y-axis. The lower branch 
equilibrium is marked with an orange cycle, the upper one 
with a black cycle, (a) For Re just slightly above the bifur- 
cation point, the basin of attraction of the upper branch is a 
small region with a bubble-like shape. Its boundary is formed 
by the stable manifold of the lower branch, (b) For Re — 188, 
the overall structure has not changed a lot, though the at- 
tractor is no longer a fixed point but a chaotic orbit, (c) At 
Re — 188.51 a boundary crisis occurs - the chaotic attractor 
becomes "leaky" and all trajectories eventually relaminarize. 
This bubble-like shape remains visible as an envelope of the 
initial conditions for longer living states, though. 
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The division by the root scales the x-axis in a way that 
the distance between LB and UB in the plot corresponds 
to the L2-norm of the difference between the two fixed 
points. 

The basins of attraction are obtained by studying the 
lifetimes, defined as the time it takes for the flow to sat- 
isfy for the first time the condition E c f(u(t) — ui am ) < 
e = 10 -8 . Initially, when the basin of attraction is closed 
and trajectories are attracted to a stable fixed point that 
is not the laminar profile, there is a solid block of long, 
actually infinite, lifetimes. The stable fixed point lies 
somewhere in the center of the blob and the unstable 
saddle lies on its boundary. 

A notable feature of these plots is the bubble shape of 
the basins of attraction of the stable fixed points. The 
unstable manifold of the unstable fixed point that loops 
around the stable fixed point folds back to approach the 
other branch of the stable manifold. This behaviour was 
identified by Lebovitz [36 , 37 in a model flow, and as the 
present figures shows it also is relevant to full represen- 
tations of the flow. In view of the bifurcation structures 
described in [29 one can expect that this bubble shape 
is rather generic. 

As Reynolds number increases, the stability proper- 
ties of LB do not change up to at least Re = 400. On 
the other side, UB is only an attractor for a very short 
range and undergoes a Hopf bifurcation at Ren = 166.05, 
which has already been noted by Clever and Busse [28] . 
The bifurcation is supercritical, as is revealed by the 
emergence of a limit cycle whose amplitude grows as 
VRe-Re H . The emerging stable periodic orbit under- 
goes a period doubling at Re = 178.9 which, for this 
domain size, results in an incomplete period doubling 
cascade. 

While the way to represent fixed points in a bifurcation 
diagram is straightforward, it is less so for periodic orbits 
or chaotic trajectories. To visualize the attracting state 
once it gets more complicated we introduce a reduced 
representation of the flow by focussing on the maxima 
in energy: we calculate a typical trajectory originating 
near UB (and hence not lying in the basin of attraction 
of the laminar state), omit an initial transient and all 
intermediate times except for the maxima. We then keep 
values of the maxima in ^E c f to represent the flow. The 
method is illustrated in figure |2j where a short section of 
a typical chaotic trajectory for Re = 187.8 is shown. The 
maxima are marked with symbols corresponding to the 
symbolic dynamics that will be introduced in section [Till 

The bifurcation cascade is summarized in figure [3j 
where the square root of the cross-flow energy E c f is 
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FIG. 2. Example of a typical chaotic trajectory at Re — 
187.8. For the reduction to a map, we only keep the energies 
at the maxima. The different symbols assigned to the maxima 
follow from the symbolic dynamics that will be introduced in 
section |Hl| squares represent the symbol "0" and circles the 
symbol " 1" . 



plotted. The attractor is plotted with a dot for every 
maximum of E c f, LB is a dashed line and UB is a dot- 
ted line once it becomes unstable. In this representation 
the emergence of a limit cycle at Ren with an amplitude 
that grows as ^jRe — Ren is clearly visible. It is fol- 
lowed by a period doubling cascade resulting in a chaotic 
attractor. We find a stable period three window from 
Re = 183.8 to Re = 185.0, followed by a new period 
doubling cascade of the period three orbit. The basin of 
attraction is still confined to a small region in state space, 
as is evident from figure [TJb), which shows the basin for 
Re = 188. Compared to the above picture at Re = 164, 
the general shape of the basin has not changed. Near 
the crisis bifurcation perhaps also the homoclinic tangles 
described in [38] form. The transition in the lifetimes is 
in agreement with the phenomenology observed also in 
the two-dimensional map described in [39 . 

For Re > Re c = 188.51, there are no more dots in 
figure [3] all trajectories eventually relaminarize. Our in- 
terpretation is that at Re c the chaotic attractor touches 
the fixed point LB (or, equivalently, its stable manifold) 
and a boundary crisis [40] occurs and the attractor be- 
comes "leaky". Even though the crisis is not visible in 
figure [3] due to the chosen representation, it will become 
visible later in another view of state space, see figures 
11 below. We show a phase space plot of Re = 195 in 
figure [TJc) . The bubble shape that formed the basin of 
attraction of the chaotic attractor before the crisis is still 
visible. But the lifetimes of the initial conditions from 
that region now are no longer infinite but vary drastically 
between almost direct decay and very long survival. The 
distribution of the lifetimes (not shown here) turns out 
to be exponential - a signature of a chaotic saddle left 
behind by the boundary crisis [25 . 



III. SYMBOLIC DYNAMICS 

In this and the next section we focus on the situation 
just before the crisis, namely at Re = 187.8, where the at- 
tractor is chaotic. We now use the reduction of the flow 
to the location of its maxima to introduce an approxi- 
mately one-dimensional map of the flow. We denote the 
value of the i-th maximum of \/E c f along this trajectory 
by Xi and plot Xi+\ vs. X{. This then gives a reduced, 
approximate map xi — > x i+ i = f(xi). Figure [4] shows a 
plot of the map obtained from a trajectory with 10798 
maxima. 

In this representation, the attractor looks very thin 
and almost one-dimensional, though a fractal structure 
can be seen perpendicular to it. In this sense it is qual- 
itatively similar to fractal attractors in well-known low- 
dimensional systems like the Henon map with standard 
parameters. The attractor consists of several branches, 
all of which show a distinctive maximum between xi = 
0.14 and X{ = 0.16, but the location of the maxima do 
not coincide. 

Given this representation of the attractor, we introduce 
a binary symbolic dynamics in the usual way: on each 
branch, a point is labeled as or 1 if it lies to the left 
or right of the maximum, respectively. The symbolic 
dynamics is shown in figure [5] for the same snapshot of the 
trajectory as in figure [2] We find that a partition cannot 
be based on a single value of Xi and that the branch on 
which it comes to lie also has to be taken into account: xq 
is a little lower than xg but assigned the symbol " 1" while 
xg is assigned "0" because of its different location relative 
to the maximum of the corresponding branch. The two 
rightmost squares in figure [5jb) lie on the branch of the 
rightmost maximum while the circle just above them lies 
on the branch of a maximum further left. 



IV. PERIODIC ORBITS 

In this section, we will describe the periodic orbits we 
found at Re = 187.8. We obtain initial guesses by scan- 
ning the map / (figure [4| for close returns along a trajec- 
tory and then use channelflows excellent Newton-Krylov- 
Hookstep algorithm [41] for finding periodic orbits. Since 
the attractor for this Reynolds number is chaotic, we ex- 
pect all orbits to be linearly unstable. For reasons of 
computational stability we restrict the search to orbits 
of symbolic period not exceeding five. 

With this method we found seven different orbits - 
they are listed, along with their properties in table [I] 
The first column contains the symbolic name of the or- 
bit, the second one the period of the orbit in time-units 
of Navier- Stokes. The Reynolds number where the or- 
bits bifurcate, either pitchfork or saddle node bifurca- 
tions, are given in the fourth column. The orbit 1 is 
created in a Hopf-bifurcation of UB, the orbits 01 and 
0111 result from period-doublings of the former one. The 
pairs of the two period-3 and period-5 orbits are the re- 
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FIG. 3. The square root of the kinetic energy yj E c f for the upper- and lower-branch solutions and the attractor coexisting 
with the laminar state. To visualize periodic orbits and chaotic attractors, we plot every maximum of E c f with a dot. The 
Nagata-Clever-Busse-state is born in a saddle-node bifurcation at Re = 163.8, we plot the lower branch with a dashed line and 
the upper branch with a line as long as it is stable and with a dotted line afterwards. The upper branch loses stability in a 
Hopf bifurcation at Re = 166.05, the emergence of a limit cycle is nicely visible. It undergoes a period doubling cascade. After 
the boundary crisis at Re = 188.51, no second attractor exists beneath the laminar state. 
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-1.70 • 10" 5 
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-0.0188 



TABLE I. Properties of the periodic orbits. The first column contains the symbolic sequence of the orbit, the second column 
the period in time-units of Navier-Stokes. The fourth column gives the Reynolds numbers where the orbits are created. The 
last four columns give the eigenvalues with the largest magnitudes. All orbits have one and only one unstable eigenvalue. 



suit of saddle-node bifurcations. To better illustrate the 
bifurcation sequence, we show the bifurcation diagram 
again in figure [6] with the bifurcation points indicated 
by dashed lines. The solid line at Re = 187.8 indicates 
the Reynolds number where the analysis of the symbolic 
dynamics takes place. 

The bifurcation sequence of the orbits is almost the 
same as the universal sequence described by Metropolis, 
Stein and Stein [42], with the only exception being that 
the period-doubling of 01 that creates 0111 takes place 
slightly after the saddle-node bifurcations in which the 
orbits 001 and 011 appear. We have not found the orbits 



0011, 0001, 00111, 00011, 00001 and 00101 - which is in 
accordance with the universal sequence for an incomplete 
cascade. 

The eigenvalues A^ of the monodromy matrix of the pe- 
riodic orbits are calculated using Arnoldi iteration [41]. 
From the eigenvalues A^ one can extract Lyapunov expo- 
nents A^ = In \Ai\/T. The four eigenvalues of largest mag- 
nitude and the corresponding nontrivial Lyapunov expo- 
nents are given in the last six columns of table [I] All of 
the periodic orbits have exactly one unstable eigenvalue 
Ai and two eigenvalues A2,3 that correspond to transla- 
tion along x and time that are theoretically strictly equal 
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FIG. 4. The chaotic attractor at Re — 187.8. It is obtained 
by calculating for a trajectory that contains 10799 maxima 
of E c f. The values of the maxima of y^Ecf are denoted Xi 
for the i-th maximum and the map is obtained by plotting 
Xi+i vs. Xi. It consists of several branches that lie close 
to each other and, for the purposes of labelling trajectories, 
can be collapsed into a one-dimensional uni-modal map. A 
binary symbolic dynamics can be introduced according to the 
relative location of the points: points that lie on a branch left 
or right of the maximum are labeled or 1, respectively. 

to 1 (translations along z are excluded by the enforced 
shift- and-reflect symmetry). The deviation of A2,3 from 
1 can be used as a measure of the accuracy of our nu- 
merics. A4 is the least stable eigenvalue of the orbit and 
measures the contraction in the directions perpendicular 
to the attractor. The positive first Lyapunov exponent 
varies between 0.0094 and 0.176, where the largest value 
comes from the shortest orbit 1. All orbits have exactly 
one unstable direction and their weakest contracting Lya- 
punov exponent A4 is of the same order, this supports the 
conjecture that the system can be approximated by a one 
dimensional map. 

The upper branch becomes unstable via a Hopf bifur- 
cation and at Re = 187.8 the complex conjugate pair of 
eigenvalues is still the only unstable ones. The values are 
Ai } 2 = 0.023 ±0.128 and the imaginary part corresponds 
to a period of T = 27r/^(A) = 48.93 which is not too far 
from the period of the orbit 1. 

The orbit is missing in table[lj An estimate of what its 
properties could be can be calculated from the 1 and the 
longer orbits in the following way: let no the number of 
zeros in the symbolic sequence and ri\ the number of ones. 
Then an estimate for the time-period of is obtained 
from an orbit composed of no and rt\ symbols "0" nd 
"1", respectively, can be obtained as follows: Since the 



FIG. 5. Symbolic dynamics for the same trajectory as in 
figure [2] squares are and circles 1. The rightmost square is 
left of the maximum and has to be labeled as despite the 
fact that it is further right than the leftmost circle. 
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0.0115 


-0.0217 


001 


53.44 


0.0052 


-0.0195 


0TT 
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-0.0247 


0111 
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0.0082 


-0.0264 


01101 


50.61 


0.0100 


-0.0229 


01111 


49.77 


0.0092 


-0.0288 



TABLE II. Estimates of some properties of the missing orbit 

0, calculated from the longer orbits. See text for details. 

properties of the long orbit are similar to a superposition 
of no and rt\ shorter segments, we can estimate the time 
from 

« (T - mT(I))/no (5) 

, and the Lyapunov exponent from 

Ag st «(AT-n 1 A 1 T 1 )/n T (6) 

they are given in table [TTJ All estimates indicate that the 
time-period of would be slightly shorter than that of 

1. While estimations for the first Lyapunov exponent 
Ai fluctuate by more than a factor of two, reflecting the 
fluctuations in the stability of the orbits themselves, the 
least stable exponent fluctuates less; it is slightly more 
contracting than those of the other orbits. 

Cobweb-plots of all orbits are presented in in figure [7[ 
While most of the orbits concentrate in the upper right 
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FIG. 6. Critical Reynolds numbers for the periodic orbit bifurcations, indicated with dashed lines. The straight line at 
Re = 187.8 corresponds to the Reynolds number where the analysis takes place. 



corner of the map, only 001 stretches out to the lowest 
and highest values. Since 001 is the only orbit with a 
sequence 00, this suggests that longer orbits containing 
00 would also stretch out to those more extreme regions 
of the attractor. 

In figures [8] and [9] we present orbits in a plane spanned 
by the cross-flow energy E c f and the friction factor C/, 
respectively. Figure [8] shows nicely the composition of 
the longer orbits as repetitions of the shorter ones. For 
example, the orbit 011 is very similar to 0111 and 01111, 
except that the latter ones have one (two) additional 
peak(s) that looks quite similar to 1. Note also how the 
orbits explore the ranges covered by a turbulent trajec- 
tory, lending support to the idea that periodic orbits can 
be used to represent all possible motions and that aver- 
ages could be calculated from the periodic orbits. 

We show some averaged quantities of the orbits in fig- 
ure 10 We use small letters for the fluctuating com- 
ponents and Uo for the base flow. The left plot shows 
(u) + Uo, the middle one (uu) and the right one (vv). 
In all three plots, the mean profiles for nearly all orbits 
impossible to distinguish among each other. Within the 
group of orbits of period less than six, the exception is 
001, which is marginally offset to the right in (b) and (c). 
We indicate the averages obtained from a chaotic trajec- 
tory with 5000 time units by crosses - they can hardly 
be distinguished from the averages of the periodic orbits. 

The motion of the walls injects energy into the fluid 
motion which in turn is dissipated by viscosity. The en- 
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FIG. 7. The periodic orbits, presented in cobweb-plots. Only 
the orbit 001 visits the extremer position on the attractor. 



ergy balance can be written as 



E = I -D, 
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FIG. 8. The root of the cross- flow energy yE c f for all 
periodic orbits. The composition of the longer orbits from 
parts of the shorter ones is nicely visible. 



s» 0.0 




0.01-0.005 

(uv) 



FIG. 10. Averaged quantities of the periodic orbits (lines) 
and a chaotic trajectory (crosses) as a function of y. (a) Aver- 
aged downstream velocity (b) Fluctuations in the downstream 
component (uu) (c) Average of (uv). 



where 



E 



D 



0.012 - 

0.010 

0.008 
0.012 

0.010 

0.008 
0.012 

0.010 

0.008 
0.012 

0.010 F 

0.008 

0.012 

0.010 

0.008 
0.012 r 

0.010 

0.008 
0.012 

0.010 

0.008 







50 



100 



150 



200 



250 



001 



011 



0111 



101101 



101111 



300 



350 



FIG. 9. The friction factor Cf for all periodic orbits. While 
Cf has maxima at roughly the same points as E c f, the sym- 
bolic dynamics can not be recognized. 



2aJ a 



( du 




du 




[dy 


^ 

y=i 


dy 


y=-i) 



In a plot of / vs. D, fixed points must lie on the 
diagonal I = D. For periodic orbits and also for the 
non-periodic chaotic motions, the center, defined as time- 
averages of / and D, will also lie on the diagonal. In fig- 
ure [TT] we present all periodic orbits and a chaotic trajec- 
tory. The points corresponding to the maxima of E c f are 
marked with circles for symbol and squares for 1 , the di- 
amond shaped symbol in the lower left corner marks LB. 
There are three points noteworthy in this figure: first, 
the symbols are clearly separated and a partition could 
be based on the location of a maximum in that repre- 
sentation. Second, the periodic orbits capture all of the 
qualitative features and loops of the chaotic trajectory. 
And third, one sees that the chaotic trajectory in the 
background is about to touch LB, the event that will 
trigger the crisis. 



V. FINAL REMARKS 

The detection of periodic orbits in this system was cer- 
tainly assisted by the fact that the flows studied here 
were confined to a small domain, restricted by a dis- 
crete shift and reflect symmetry, and also analyzed at 
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a) n „ a low Reynolds number. In this sense, it is similar to 

the study of fluctuations in the Kuramoto-Shivashinsky 
equation [32 . Nevertheless, the reduction to a symbolic 
dynamics compatible with a unimodal one-dimensional 
map and the association of the orbits to the univer- 
sal Metropolis -Stein- Stein sequence shows that also this 
high-dimensional system follows the universal dynamics 
identified in low-dimensional systems. This raises the 
hope that bifurcation analyses and periodic orbits can 
also be identified and used in other systems, such as pipe 
flows or plane Poiseuille flow. The ultimate goal, how- 
ever, is to analyze not only chaotic but also some turbu- 
lent flows in terms of periodic orbits, as suggested e.g. in 

El- 



"1.5 1.6 1.7 1.8 1.9 2.0 2.1 2.2 2.3 2.4 
1 




ACKNOWLEDGEMENTS 



I 




FIG. 11. The periodic orbits in the input-dissipation plane. 
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